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^Maps  of  the  monthly  median  of  refractive  effects  are  useil  by  precision 
navigation  and  radar  systems  to  lorrect  for  measurement  errors  caused  by 
the  ionosphere.  Calculated  off  line  using  available  climatologies  of  the  per- 
tinent ionospheric  parameters,  they  are  used  in  real  lime  systems  to  removr 
70  to  75  percent  of  the  monthly  rms  error.  When  greater  nrecision  is  rerpiire( 
these  maps  may  he  adapted  by  any  of  several  local  measurements  of  current 
ionospheric  conditions.  The  residual  refractive  error  is  dependent  on  the 
density  of  these  local  measurements  in  space  and  time  and  on  their  ability  to 
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Adaptive  Mapping  af  Mid-Lotitude  Ionosphere 


1.  tvvMt)tn»n'u» 


As  navig.-ition  and  radar  systems  pvolvp,  :h<*  requirements  for  precision  ranne 
measurements  become  more  sfrincent.  Kventuall  • real  time  algorithms  must  be 
developed  to  provide  increased  accuracy  in  corrections  for  refractive  effects  arising 
from  the  ionosphere.  Helatively  simple  look-up  tables  of  the  expected  median  cor- 
rection can  he  derived  from  climatologies  of  the  principal  ionospheric  parameters. 
Oependent  on  onlv  smoothed  solar  conditions,  these  tables  can  be  prepared  several 
months  in  advance,  and  when  used,  will  typically  remove  SO  to  75  percent  of  the 
monthly  rms  refractive  error.  Kor  more  precise  estimates  of  refractive  effects, 
maps  can  be  developed  from  climatology  for  the  enti  re  system  field  of  view.  Stored 
in  the  system  processor,  these  maps  may  be  adapted  by  pertrnent  local  measure- 
ments to  provide  a real  time  correction  that  will  represent  the  response  of  the 
ionosphere  to  current  solar-terrestrial  conditions  and  even  reflect  local  features, 
in  some  cases  (I)uLong  et  al>.  With  the  adaptive  techniques,  the  residual  refractive 
error  is  dependent  on  the  information  within  the  locaL  measurement  and  the  density 
of  measurements  in  space  and  time.  When  the  highest  precision  is  required,  tech- 
niques have  been  developed  for  sophisticated  systems  to  directly  measure  the 

(Heceived  for  publication  29  .luly  1977) 

t.  Dul.ong.  I). I).,  and  Allen.  H.S.  (1976)  Specification  of  Navigation  and  Radar 
Krrors  Caused  by  the  Ionosphere.  ('fisi’AR  .Symposium  r*roceedings.  l"he 
(^.eophysical  Use  of  Satellite  Vleacon  Observations,  lloston  I’niversity. 


needed^refrrictive  correction  for  the  lononphere  alone  the  path  to  the  target  (Alien 
et  al>.  This  paper  will  examine  adaptive  techniques,  using  measurements  fmm 
the  Navy  Navigation  Satellite  System  (NNSSl  to  illustrate  some  of  the  promise  and 
some  of  the  problems. 

2.  THK  RKKR  XtmVK  CORRKtTIUN 

At  I HF  and  higher  radio  frequencies  (f>  the  phvsics  for  wave  refractive  effects 
can  be  greatly  simplified;  the  local  electron  densitv  (n»  is  the  onlv  variable  that 
must  be  considered.  Then  the  ray  path  is  approximated  by  a straight  line  with  a 
first  order  correction  for  change  in  apparent  range  (All!  which  is  directly  propor- 
tional to  the  electron  content  integrated  along  the  slant  p;ith  (s). 
s 

AK  = k/f^  I nds  . ,,, 

0 

3 

Millman  has  shown  that  from  this  a correction  can  be  made  for  the  first  order 
change  in  elevation  angle.  For  targets  in  the  upper  ionospliere  and  at  low  elevation 
angles  the  change  in  the  elevation  angle  (AK)  is 


AE  = (cotE^)  Alt/H  . 

where  F.o  is  the  initial  elevation  angle.  The  principle  correction.  Alt.  may  be 
modeled  from  existing  climatologies  of  ionospheric  pa rameters  (t.lewellvn  and  Hent*.'* 
Such  a modeling  shows  that  the  median  of  the  range  correction  induced  bv  the  ionos- 
phere has  greatest  values  just  after  midday,  at  the  equinoxes,  for  conditions  at  the 
maximum  of  the  solar  cycle.  This  correction,  in  meters,  is  illustrated  in  Figure  1 
for  a representative  radar,  using  425  MMr  along  a slant  path  to  a target  with  eleva- 
tion angle  S deg  and  altitude  1000  km.  These  curves  represent  the  twelve  month 
running  mean  sunspot  number  CRp  at  the  maxima  of  the  last  20  cvcies  of  solar 
activity  when  the  peak  TT^  of  a cycle  could  range  from  SO  for  a verv  quiet  to  200  for 
a very  disturbed  solar  maximum.  The  fractiles  of  TF  are  the  lower  quartile.  med- 
ian.  upper  qua rtile.  and  upper  decile  (70.  100.  1.10.  and  ISO  respectively!.  If  the 

2.  Allen.  R.S. . Duj.ong.  D.  I). . Grossi.  M.l). . and  Kat*.  A.  II.  (1076)  Ionospheric 

Error  Correction  in  Precision  Radar  Systems  by  Adaptive  n'roTilng 

of  The  Fropagation  Medium.  AGARD  TonfereAce  Proceedings.  Pcnpaftnlfi^n 
Cl^OO^'on^  of  Navigation  and  Positioning  Systems.  .Ankara.  Turkey.  AGARI) 

3.  Millman.  G.  H. , and  Reinsmith.  G.  M.  (1074)  An  Analysis  of  the  Incoherent 

Scatter- Faraday  Rotation  Technique  for  lonospnenc  Propacation  Frror — 
f orrection.  (.eneral  Electric  Co.  Tf?^  1174 1'MTlg. 

4.  l.lewellyn.  S.K..  and  Rent,  R,  H.  (1071)  Documentation  and  Description  of  the 

Rent  Ionospheric  Model.  A FGRl .-TR-73-06R7, — '' 
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mapping  of  a correction  is  made  for  the  entire  field  of  view  of  a avatem  uiiing  param* 

eters  in  the  clin.at'^logy  that  represent  actual  me<iian  ionospheric  conditiona.  when 

it  is  used  to  correct  real  time  range  measurements,  tlie  remaining  monthly  rma 

error  (6K  ) is  al>out  ’JO  to  2!>  percent  of  daytime  and  al<out  10  to  IS  percent  of 
*’1  S 

nighttime  median  v.ilues  (Itul.ong).  ‘ For  the  r.idar  of  l-  igure  I.  an  eaperted  Alt  at 

the  next  solar  maximum  near  the  equinoxes  is  300  m,  resulting  in  a Alt  of  about 

m 

70  meters.  Systems  needing  greater  precision  mav  further  re«luce  refractive  effects 
by  using  some  local  ionospheric  measurement  to  adapt  this  median  prediction,  in  a 
manner  similar  to  that  described  bv  l>ul.ong  et  al.  * 
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l-  igure  I.  Monthly  .Mean  tlnngr  t'orrortion  at  Middav 
for  a 425  Mil*  K idar  Situ.ated  Near  50*\  l.atitude. 
Target  is  at  5®  elevation  angle,  100  km  altitude. 

Kach  curve,  parametric  on  12  month  running  mean 
sunspot  number,  represents  tlie  expected  range 
correction  for  solar  maximum  conditions.  \t  an 
average  solar  maximum.  H = * 100,  the  expecteil 
correction  is  near  300  m it  the  equinoxes 
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The  residual  error  when  using  an  adapted  model,  AH_^,  depends  on  the  apace 
and  time  coherence  of  day-to-day  changes,  including  local  features  which  occur 
within  the  r:idar  field  of  view.  To  examine  this  error,  actual  TKC  data  was  ob- 
tained from  observations  of  N'NSS  satellites  for  four  days  in  December  1075,  taken 
at  l.indau,  I-  H(7,  and  (Ira?,  .Austria,  and  for  three  davs  in  October  1074,  from 
Bozeman  and  Boulder  in  Colorado  (1  eitinger  and  Hartmann).*’  These  represent 

5.  Dul.ong,  D.  I).  (1077)  Ueduction  of  the  Tncertainty  of  lladar  Mange  Correction. 

.A FCil .-TH-77  (in  press'. 

6.  I. eitinger,  II.,  and  Hartmann,  Ci.K.  (1076)  Time  and  l.atitude  Dependence  of 

Ionospheric  Klectron  Content  from  the  Combination  of  NN^S  ann  AT^-(i  Hatn, 
CO.''1’aI(  .Symposium  Proceedings,  the  Geophysical  Use  of  ?(atellite  Deacon 
Observ.itions,  Boston  I.'nivcrsity. 


quiet  and  disturbed  periods  of  K^'^^Kuetic  activity,  respectively.  The  observations 
were  normalized  to  simultaneous  TKC  data  reduced  from  single  station  observa- 
tions of  the  beacons  of  the  ATS-6  satellite.  Hv  converting  the  NNSS  data  to  radar 
range  correction,  the  results  can  be  used  to  adapt  maps  of  the  median  range  cor- 
rection and  the  reduction  in  refractive  error  can  he  examined. 

In  simulating  a proce<lu“e  for  adaptive  modeling,  data  from  the  first  in  a 

sequence  of  NNSS  passes  may  be  use<!  to  normalize  the  prediction  of  median  range 

correction  to  conform  to  the  actual  conditions  at  that  time.  The  scaiing  used  along 

the  latitudinal  interwils  for  this  pass  would  then  be  used  to  adapt  the  median  along 

the  path  at  the  time  of  the  next  NNSS  pass.  I'he  first  pass  may  be  considered  a 

representation  of  a local  ionospheric  measurement  and  the  subsequent  passes  may 

represent  system  measurements,  such  as  a radar  svstem  measuring  the  range  to 

a target.  The  error  Oil  . can  be  compared  to  the  error  flit  . to  determine  the 

m 

effectiveness  of  the  adaptive  modeling  procedure.  This  may  be  repeated,  using 
each  pass  to  adapt  the  model  and  comparing  the  predictions  of  the  adapted  model  at 
the  time  of  subsequent  NNSS  passes  with  the  actual  observations. 

Some  of  the  effects  of  using  this  procedure  to  specifv  refractive  errors  in  day- 
time can  be  seen  in  f igure  2 for  the  quiet  period  in  December  1?>7S.  The  upper  set 
of  curves  illustrate,  on  a logarithmic  scale,  the  range  correction  that  the  system 
would  have  made  along  the  track  of  an  NNSS  pass  at  1240  I.T  on  12  December  as 
determined  from: 

(D  the  NNSS  observ-ations  themselves,  dotted  line. 

(2>  the  original  median  predicted  by  climatology,  dashed  line. 

(31  the  me<tian  scaled  by  the  prior  NNS.S  pass  at  1 120  I.T.  solid  line. 

The  lower  set  of  curve  are  on  a linear  scale  and  represent  the  absolute  value  of 
the  residual  range  error  occurring: 

(11  using  the  median  from  climatology.  dashed  line. 

(21  using  the  adapted  median.  6H  . solid  line. 

The  gradients  within  the  predicted  median  do  not  match  those  of  the  NNSS  observa- 
tions, creating  large  errors  in  some  regions,  but  these  are  reduced  significantly 
by  the  adapted  prediction.  Note,  however,  that  using  all  of  the  information  along 
the  1 120  pass  predicts  local  features  which  are  either  absent  or  displaced  at  the 
time  of  the  next  pass.  Similar  effects,  observed  during  the  (ictober  1974  period 
are  shown  in  Figure  3 for  the  da>-time  hours  of  12  (Ictober.  The  observed  values 
of  radar  range  correction  along  the  latitudinal  track  of  consecutive  NNSS  passes  at 
0920,  1040,  and  1550  I.T,  solid  lines,  rapidly  become  enhanced  with  respect  to  the 
predicted  monthly  median,  dashed  lines,  for  these  local  times.  If  all  of  the  informa- 
tion along  the  0920  I.T  pass  were  used  to  scale  the  median  in  order  to  better  pre- 
dict the  1040  I.T  ionosphere,  and  similarly  the  1040  pass  were  used  to  specify  the 
1550  ionosphere,  two  effects  would  occur.  First,  since  the  general  enhancement 
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persists  over  several  hours,  the  systematic  bias  of  pre<liction  would  be  reduced. 
Second,  since  ma{or  local  features  with  amplitude  of  about  10  percent  have  either 
move<i  in  space  or  have  changed  character  with  time,  there  is  an  additional  error 
that  indicates  smoothing  may  be  necessary  in  tlie  adaptive  modeling. 
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Figure  2.  Adaptive  Modeling  I'sing  Data  From  NN'SS 
Measurements.  The  first  pass  at  1120  I.T  is  used  to 
scale  the  predicted  latitude  variation  of  the  range 
correction  at  1120  l,T  and  then  the  scaled  values  are 
predicted  forward  to  1240  I.T.  The  envelope  of  the 
absolute  error  using  the  adapted  median  is  significantly 
less  than  using  just  the  monthly  median.  Note  that 
local  features  along  one  pass  do  not  match  local  features 
of  the  other  hut  that  both  have  a consistent  ’atitude 
variation  quite  different  from  that  of  the  median 
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I'igure  3.  General  Knhancement  of  the  Hange  Gorrertions 
Over  a Wide  Latitude  Mange  on  12  October  1974.  The 
systematic  bias  between  prediction  and  observation  could 
be  removed  by  adapting  the  median  prediction  with  simple 
scaling  factors.  Local  features  of  about  10  percent 
amplitude  should  be  filtered  out  if  the  adaptions  are  spaced 
an  hour  or  two  apart 


During  the  early  morning  hours  of  this  enhanced  day,  <12  October)  the  ionos- 
phere was  depressed  below  the  expected  median  in  some  regions,  Ligure  4.  By 
0400  I.T  there  was  a deep  trough  in  TKC  centered  near  49“N  Geographic  (L=41  and 
a diffuse  enhancement  in  the  plasmapause  region  near  40°N  (L=2.  7).  At  0R40  I T 
these  had  both  moved  southward  about  2°  and  the  relative  depth  of  the  trough  had 
decreased  while  the  plasmapause  enhancement  had  increased.  The  smoothed  pre- 
sence of  these  features  may  be  discerned  in  the  expected  monthly  medians  derived 
from  climatology.  Note  that  the  persistent  depression  of  observed  values  during 
this  night  was  not  carried  into  the  sunlit  hours;  by  0920  I T the  oh.cervations  were 
almost  at  the  expected  values.  This  is  consistent  willi  a previous  study  of  the 
temporal  variability  of  radar  range  error  at  a fixed  station,  in  which  it  was  shown 
that  useful  predictions  could  not  he  made  across  the  sunlit  terminators  (DuLong',' 
The  early  morning  of  12  October  1974,  was  a relatively  quiet  period  preceding 
the  magnetic  disturbances  that  affected  the  major  portion  of  the  October  sample. 
The  features  observed  during  this  period  were  also  observed  for  each  of  the  late 
afternoon  through  early  morning  periods  of  the  quiet  sample  from  December  1975. 
In  the  late  afternoons  of  each  of  the  December  days  there  is  an  enhanced  region 
near  RO^N  (1,5:3.  5)  that  is  a well  developed,  persistent  feature,  shown  in  the  upper 
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set  of  curves  in  KiKurc  5.  Hy  backtrarkinj;  from  2000  LT  in  the  central  set  of 
curves,  where  the  trough  in  TKC  near  55“N  (I.s5l  can  be  identified  each  day,  this 
feature  in  the  December  afternoon  can  be  identified  us  a filasmapuuse  enhancement 
of  TKC.  The  lower  set  of  curves  in  Kigure  S show  that  both  the  plasmapause 
enhancentent  and  the  trough  persist  into  the  early  morning  hours,  as  was  also 
observed  in  Figure  5 for  the  October  tlay. 
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Figure  4.  Sample  of  the  Fffect  of  the  Trough  on 
Flectron  Density.  Seen  Near  49“ N Latitude  on 
12  October  197  4,  on  the  Fxpected  Hange  Correction 
for  a 425  MMz  Kadar.  Predictions  adapted  for 
this  feature  were  valid  over  severa'l  hours  on 
this  night.  Contrast  the  general  depression  with 
respect  to  the  predicted  median  to  the  general 
enhancement  seen  the  following  morning  (Figure  3) 
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l-'imirc  .S.  Sample  of  Cniiseeiitive  Days  When  the 
Observations  of  NNSS  Showed  that  the  Position  and 
Charaoteristies  of  the  I'roiiL'h  Near  fi.‘i'’N  l.atitnde 
anti  the  Plasmapause  Pnhaneement  Near  SO^N 
I.atituiie  Could  be  Sperified  f rom  lIour-to-Hour 
and  1 lay -to -Dav,  Ad.apted  predietinns  scaled 
from  these  observations  could  reiluce  the  error 
of  rani>e  correction  bv  several  meters 


As  further  illustration  of  what  may  be  encountered,  consider  the  field  of  view 
of  a hyixithetical  radar  situated  at  80*W,  Long.,  40'’N.  Lat. , (Figure  61.  The  range 
correction  along  a latitudinal  swath  near  the  western  horizon  can  be  determined 
from  \NSS  passes  observed  at  Boulder,  Colorado,  as  shown  in  the  sample  from 
October  1974.  In  actual  system  use,  these  would  reflect  the  increase  in  slant 
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Figure  6.  Coverage  Diagram  of  a Hypothetical  Radar  That  Illustrates  the 
Need  to  Determine  How  liata  From  One  Part  of  the  Field  of  View  May  Rest 
be  Used  to  Adapt  Predictions  for  the  Remaining  Portion.  Continuous 
measurements  of  total  electron  content  derived  from  ATS-3  observations 
along  the  eastern  seaboard  can  be  compared  to  NNSS  observations  made 
at  mid-continent 
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thickness  of  the  ionosphere  as  elevation  angle  decreaseH.  but  here  it  is  convenient 
to  stay  with  equivalent  vertical  values  to  be  consistent  with  the  reduced  NNSS  data. 
Now  the  question  that  arises  is  whether  adapted  features  can  be  mapped  across  the 
entire  field  of  view.  Using  single  station  TKU  measurements  obtained  using  the 
\ IIK  beacon  on  the  ATS-.'l  geostationerv  satellite  at  Narssarssuaq,  Greenland; 

Goose  Hay.  l abrador:  Hamilton,  Massachusetts:  and  Kennedy  Space  Flight  Center, 
Florida,  it  is  possible  to  construct  a similar  latitudinal  swath  near  the  eastern 
horizon  of  the  radar.  1 he  features  of  this  swath  can  then  tie  compared  with  the 
NN.SS  data  near  the  western  horizon.  A comparison  of  the  data  in  I.  value  and  local 
time  for  3 sequential  nighttime  passes  on  12  October  1974.  a quiet  period,  and 
3 sequential  daytime  passes  on  13  October  1974,  a magnetically  disturbed  period, 
are  displayed  in  Figures  7 and  8,  respectively.  On  the  quiet  night,  the  sequence  of 
development  of  the  trough  in  local  time  near  the  eastern  horizon  and  its  character 
and  position  in  1,  value  are  reproduced  near  the  western  horizon.  It  is  then  possible 
to  conclude  that  it  would  have  been  useful  to  map  it  over  the  whole  field  of  view.  In 
contrast,  during  the  disturbed  day  a large  enhancement  of  the  ionosphere  developed 
over  the  eastern  portion  of  the  radar  field  of  view,  hut  no  similar  feature  occurred 
in  the  western  portion.  At  I 100  I T on  13  October.  TKC  at  Goose  H.ay  was  already 
higher  than  at  adjacent  latitudes.  Hy  1200  I T TKC  at  Hamilton  and  Goose  Hay 
increased  slightly,  but  at  Narssarssuaq  it  increased  significantly.  At  HiOO  l.T  TKC 
at  Narssarssuaq  had  decreased  to  the  value  observed  at  1 100  l.T;  at  Hamilton  there 
was  a further  decrease,  while  at  Goose  Hay  there  was  an  extreme  enhancement.  On 
the  western  horizon  of  the  radar  there  was  an  increase  near  1200  l.T  in  the  northern 
coverage  of  the  NNSS  passes  that  agrees  with  the  enhancement  seen  at  that  local 
time  at  Narssarssuaq.  hut  nothing  in  the  western  data  is  comparable  to  the  major 
enhancement  observed  at  Goose  Hay.  During  such  greatly  disturbed  periods  adaptive 
modeling  would  be  of  most  use  to  a system  but  this  example  suggests  that  extensive 
mapping  of  local  features  in  either  latitude  or  longitude  should  not  be  undertaken 
without  sufficient  knowledge  of  their  spatial  and  temporal  characteristics. 
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RAOAR  RANGE  CORRECTION.  RADAR  RANGE  CORRECTION 


Kiijure  7.  Dcvplopnient  of  tlio  Trough  Alont;  the 
Kastorn  Seaboard,  Determined  from  ATS -.3 
Observation  at  Narssarssiiaq,  (kiose  Hay, 

Hamilton,  ami  Kennedy  (N,  (1,  H,  and  K) 

Compared  With  the  Development  Over  the  Mill- 
Continent  as  Determined  From  N’NSS  Observations. 
For  this  quiet  ni^ht,  larije  features  may  be  mapped 
m loral  time  and  macnetic  coordinates  over  the 
entire  radar  field  of  view 
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Figure  8,  Same  Coverages  as  Figure  7.  For  this 
daytime  sample  there  was  a very  large  enhancement 
of  the  ionosphere  at  high  latitudes  that  was  confined 
to  the  eastern  segment  of  the  field  of  view.  During 
such  periods,  features  could  not  be  adapted  across 
the  entire  field  of  view 
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Actual  observations  made  with  the  Navy  Navication  Satellite  Systems  (NNSS) 
satellites  have  been  used  to  test  a technique  which  would  adapt  maps  predicting 
ionospheric  refractive  effects  over  the  field  of  radar  or  navigation  systems.  In 
this  adaptive  technique  each  NNSS  pass  is  used  to  scale  local  features  into  the 
monthly  median  maps.  Kach  adapted  map  is  then  used  to  predict  the  local  behavior 
of  the  refraction  correction  and  the  results  compared  with  ensuing  NNSS  passes. 

Systematic  changes  which  represent  the  daily  response  of  the  local  ionosphere 
to  solar-terrestrial  conditions,  such  as  the  general  enhancement  seen  over  the 
Boulder  region  on  12  October  1074,  are  usefully  predicte<i  by  this  technique  but 
with  the  restriction  that  a prediction  should  not  he  carried  across  either  sunlit 
terminator  as  this  could  increase  error. 

During  quiet  geomagnetic  periods,  those  local  features  with  large  geophysical 
scale  size  and  a consistent  identification  in  climatology,  such  as  the  nighttime 
trough  and  the  plasmapause  enhancement,  seen  so  persistently  in  the  10-13  December 
TEC  observations  from  Europe,  may  be  profitably  tracked  by  local  measurements 
and  indeed  may  he  usefully  extended  over  the  field  of  view.  During  the  viavdime  some 
local  features  may  be  as  large  in  amplitude  but  so  transitory  in  either  space  or  time 
that  thev  are  not  observed  on  subsequent  passes.  In  this  case  it  is  better  to  smooth 
these  unknown  features  prior  to  adapting  the  median  maps. 

.As  is  well  known,  during  ciisturbed  geomagnetic  conditions  some  local  features 
maybe  very  dynamic.  Therefore,  they  must  be  closely  specified  in  the  adaptive 
technique  during  such  periods.  Moreover,  very  large  changes  in  the  ionosphere 
can  occur  in  constrained  portions  of  the  coverage  area,  such  as  the  intense  enhance- 
ment observed  on  13  October  107  4 over  the  eastern  seaboard  of  North  .America  but 
not  at  the  mid-continent  stations.  Such  features  must  be  carefully  defined  in  space 
and  time  and  can  not  be  extended  beyond  the  region  of  observation. 
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